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Living organisms process information without any central control unit and without any ruling 
clock. We have been studying a novel computational strategy that uses a geometrically arranged 
excitable field, i.e., "field computation." As an extension of this research, in the present article we 
report the construction of a "direction detector" on an excitable field. Using a numerical simulation, 
we show that the direction of a input source signal can be detected by applying the characteristic as 
a "coincidence detector" embedded on an excitable field. In addition, we show that this direction 
detection actually works in an experiment using an excitable chemical system. These results are 
discussed in relation to the future development of "field computation." 

PACS numbers: 05.45.-a, 05.70.Ln 



I. INTRODUCTION 

Modern digital computers process information in a 
step-wise manner using a CPU equipped with a ruling 
clock. Improvements in the step-wise operation, with re- 
gard to speed and reliability, have made digital comput- 
ers indispensable in modern society. However, step-wise 
operation under the control of a CPU inevitably exhibits 
a serious bottleneck when performing "real-time opera- 
tions" on large sets of time-dependent input data. To 
break through such a bottleneck, many approaches to 
"parallel computation" have been proposed. Most previ- 
ous studies on "parallel computation" seem to have fo- 
cused on the problem of how to effectively connect multi- 
ple unit computers (see, e.g., [lj]). On the other hand, liv- 
ing organisms routinely perform "parallel computation" 
in the absence of any ruling clock. Higher order animals 
obtain the information from the external world through 
sensory systems, make decisions, and act optimally, in 
"real time" and in the proper order. Thus, an investiga- 
tion of "parallel computation" that mimics the essentials 
of computation in living organisms may be of scientific 
value. 

It is well known that characteristics of excitability 
can be found in a wide variety of biological systems, 
including nerve membranes, and that living organisms 
use excitability in cells or tissue for information process- 
ing Hi HI El Studies on information processing in ex- 
citable media may lead to a fundamental understanding 
of information processing in biological systems. Several 
studies have provided interesting examples of informa- 
tion processingwith an excitable field, such as in image 
processing 0, IS 0, HI an d logic operations 0, 0, ^1 • 
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Recently, it has also been reported that, in both a simu- 
lation 0,0,0,0,03 and an actual experiment on an 
excitable system |l7l Il8j , "diode" can be made to appear 
between opposing excitable fields by choosing the spatial 
arrangement of excitable fields in an asymmetric manner. 

By applying the unique characteristics of such diodes 
on an excitable field, logic gates with the operations 
AND, OR, and NOT can be created 0, where the ef- 
fect of the coincidence of a pair of signals with a time 
difference plays an essential role. Such logic gates nat- 
urally lead to "time-dependent" signals in the absence 
of any timing apparatus, i.e., without a "clock." As an 
extension of these logic gates, interesting examples of in- 
formation processing have been reported such as a time- 
interval detector, a comparison detector, and a subtrac- 
tion detector 0|- It nas also been reported that, with ex- 
citable fields in a comb-like arrangement, the coincidence 
or time-difference detector of a couple of input pulses 
can also be constructed in a real experiment using an ex- 
citable chemical s yst em the Belousov- Zhaboti nsky 
(BZ) reaction [20H21I]. and in computer simulations |l4j . 
The detection of a time-difference between a pair of in- 
put pulses is one of the basic functions of information 
processing in nervous systems, such as in sound localiza- 
tion in auditory systems [22, |23| . 

In this study, by extending the idea of a "coincidence 
detector," we report a simulation, and a corresponding 
actual experiment, of a "direction detector" that can de- 
tect the direction from which a chemical wave propa- 
gates, using an excitable field with a suitable geometric 
arrangement. 



II. NUMERICAL METHODS 

As an actual experimental system, in the present study 
we adopted the photosensitive BZ reaction. In this chem- 
ical system, illumination with light results in the produc- 
tion of bromide, which inhibits the oscillatory reaction. 
Thus, the excitability can be adjusted by changing the 



2 



intensity of illumination 0, fl3L 01 . The essential fea- 
tures of the photosensitive BZ system can be described by 
a two- variable Oregonator model [2^ |2|j that has been 
modified to include the photochemical pathway j23, , 
as in Eq. ^ 
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where the variables it and f correspond to the dimen- 
sionless concentrations of activator and inhibitor, respec- 
tively. D u is the diffusion coefficient of the activator. The 
parameter <j) corresponds to the light intensity, and the 
local excitability of the medium is inhibited when <fi is 
high 0. The kinetic parameters e = 0.05, q = 0.000f5, 
/ = f .0 and D u = 1.0 are taken to be constant through- 
out the simulation. The medium is excitable at <f> — 0.007 
and becomes nonexcitable as 4> increases. Numerical sim- 
ulations were carried out with Eq. using the ADI 
method (alternating direction implicit method), with dif- 
ferential calculus and the Runge-Kutta method 30] . The 
boundary condition at the edge of the frame is taken to 
be no flux (Neumann condition), while the flux between 
the excitable and inhibitory fields is taken to be free. The 
grid size is 300 x 300 points in a square lattice, and the 
time interval is taken to be At = O.OOOf . The mesh size 
is taken to be 0.2. 



III. NUMERICAL RESULTS 

Figure ^ depicts a "direction detector" constructed of 
excitable and inhibitory fields. In Fig. [J a), both the 
deep and shallow purple regions correspond to excitable 
fields (</> = 0.007 and <j> = 0.033 in Eq. QJ, while the white 
region corresponds to an inhibitory field (tp — 0.209). 

Figures Wfo) and Wi c ) show time series of the propa- 
gation of a signal wave. In Fig. ^b), a wave triggered 
at point Si (Fig. QJb) t —1 ) enters the ring-shaped ex- 
citable field through the point P\ and splits into a pair 
of rotating pulses (Fig. ^b) t =4), which causes the sig- 
nal to be transmitted to output channel X 2 (Fig. ^b) 
t =10). On the other hand, in Fig. Q^c), a wave trig- 
gered at point £2 (Fig. Wi c ) t =1) enters the excitable 
ring through point Pi and splits into a pair of rotating 
pulses (Fig. \T^c) t =6), which causes the signal to be 
transmitted to output channel X±. Thus, the transmis- 
sion of the signal onto an output channel is generated 
only if the pair of pulses collide with each other near the 
gate to the output channel. It is clear that the output 
channels X1-X3 provide the information on the original 
direction of the input signal. 
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FIG. 1: (a)Schematic representation of a "direction detector" 
on an excitable field. The shallow and deep purple regions cor- 
respond to excitable fields, while the white region corresponds 
to an inhibitory field. The field geometry is arranged so that 
a signal wave is transmitted into one of the output channels 
X\ — A3 depending on the direction that a wave triggered in 
the shallow purple area propagates. (b)An example of a prop- 
agating wave. A wave triggered at point Si causes a signal to 
be transmitted to output channel A2. (c) Another example of 
a propagating wave. A wave triggered at the point S2 causes 
a signal to be transmitted to output channel Ai . 



IV. EXPERIMENTAL METHODS AND 
RESULTS 



We examined the above idea in a real experiment. 
As an experimental model, we adopted a photosensi- 
tive version of the BZ reaction under conditions sim- 
ilar to those in previous studies 0, 0]. The water 
was purified with a Millipore-Q system. The BZ solu- 
tion we used in this experiment consisted of 6.8 x 10 _1 M 
sodium bromate (NaBrOa), 3.0 x 10 _1 M sulfuric acid 
(H 2 S0 4 ), 2.0 x lO^M malonic acid (CH 2 (COOH) 2 ), 
5.0 x 10~ 2 M sodium bromide (NaBr), and 1.7 x 10" 3 M 
ruthcnium(II)-bipyridyl (Ru(bpy)g + ). Cellulose- nitrate 
membrane filters (Advantec, A100A025A) with a pore 
size of 1/im were soaked in BZ solution for about 1 
minute. The membrane was gently wiped with filter pa- 
per to remove excess water and placed in a petri dish. 
The surface of the membrane filter was immediately cov- 
ered with silicon oil (Shin-Etsu Chemical Co.) to pre- 
vent it from drying and to protect it from the influ- 
ence of oxygen. The experiments were carried out in 
an air-conditioned room at 21 ± 1°C, at which the re- 
action medium showed no spontaneous excitation. Wc 
created a geometry composed of white, gray, and black 
regions, similar to the circuit in Fig. ^a), with a personal 
computer. Such image was illuminated on the medium 
from the below through a liquid-crystal projector (Ep- 
son, ELP-810) as shown schematically in Fig. [5] The 
liquid-crystal projector was used as the light source, and 
the spatial intensity distribution was controlled by the 
computer. The experiments were monitored from above 
with a digital video camera (Panasonic, NV-DJ100). For 
image enhancement, a blue optical filter (AsahiTechno- 
Glass, V-42) with a maximum transparency at 410 nm 
was used. The light intensity at the illuminated area was 
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determined by a light intensity meter (ASONE, LM-332). 

Figure [3] shows the experimental results in the pho- 
tosensitive BZ reaction, where the white and purple re- 
gions correspond to inhibitory and exciteble fields, re- 
spectively. In Fig. |2| the white, shallow purple, and 
deep purple regions are illuminated at the intensity of 
6.4 x 10 4 lx, 6.6 x 10 3 lx, and 2.0 x 10 3 lx, respectively. 
Chemical waves are colored orange. This detector is 
made as simple as possible, so it has only one output 
channel (A 2 ). In Fig. Ola), a chemical wave triggered at 
the point Si propagates into the ring-shaped excitable 
field through point Pi , which is the nearest to S\ on the 
excitable ring. The wave then splits into two pulses that 
undergo clockwise and counterclockwise rotation, respec- 
tively. The two pulses collide with each other and disap- 
pear at point Ri, which is opposite P\, and a signal is 
transmitted to output channel X 2 (Fig. Ota) 223s). In 
a similar way, in Fig. |3Jb), a chemical wave triggered at 
S2 enters the ring through P 2 • The wave then splits into 
two pulses that undergo clockwise and counterclockwise 
rotation and collide with each other at point i? 2 , which 
is opposite P 2 (Fig. Elb) 228s). No signal is transmitted 
to output channel A 2 , even though the clockwise pulse 
passes near the gap between X 2 and the excitable ring. 
Thus, the essential characteristics of the "direction de- 
tector" in the numerical simulation were realized in the 
experiment. 




FIG. 3: Experimental results with a "direction detector" us- 
ing the photosensitive BZ reaction. This detector is the sim- 
plest version with only one output channel (X2). In order to 
show the result as an illustration manner, quasi-color is used 
in the figure. (a)A signal wave is transmitted to the output 
channel X2 when a wave propagates from above. (b)No signal 
wave is transmitted to output channel X2 when a wave does 
not propagate from above. 
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FIG. 2: Scheme of the experimental setup. 



V. DISCUSSION AND CONCLUSION 

We have shown that, in both a simulation and an ac- 
tual experiment on an excitable system, a "direction de- 
tector" can be created by suitbly arranging the spatial 
geometry of excitable fields. Although we adopted the 
Oregonator (Eq. ^) for an excitable field in the numeri- 
cal simulation in our present study, the essential feature 
of wave propagation on an excitable media can also be de- 
scribed by the FitzHugh-Nagumo equation [3lH32| . since 



both are reaction-diffusion equations. In our previous 
studies, we have shown various kinds of "field computa- 
tions" in numerical simulations using FitzHugh-Nagumo- 
type equation [TiLIT^]. 

In the "direction detector" as Fig. ^a), the location 
of the azimuth is converted into a timing difference be- 
tween gates A and B. Next, through the AND gates 
(logic product) at the output Xk (k = 1,2,3), this tim- 
ing difference is detected as the position of the collision. 
Thus, part of the mathematical structure of the "direc- 
tion detector" can be easily described by using a function 
as in Eq. [21 



X[o,a]K x > - 1 1 (xe[0,a\). 



(2) 



Let us suppose in Fig. ^a) , after a wave is trig- 
gered at a certain time and in a certain direction 9, two 
pulses are triggered at the points A and B at time t\ 
and i 2 , respectively. That is, the location of the az- 
imuth is converted into information in the two pulses 
which travel along the ring-shaped excitable field with a 
time difference. If we associate the arrival of one pulse 
with an input of logical true, then the input to gate Xk 
(k = 1,2,3) by the counterclockwise pulse is described 
by X[o.a] ( c (t ~ ti) ~ irdk/n), where d ~ AO = BO, c is 
the propagation speed of the pulses on the excitable ring, 
and n is the number of the segments divided by the out- 
put channels Xk, i.e., the number of output channels +1 
(n — 4 in Fig. H (a)), a is a positive constant and the 
interval [0, a] indicates the permissible range of the time 
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difference: the excitability of the field and the shape of 
output gates are reflected in this parameter. In the 
same way, the input to gate Xk by the clockwise pulse 
is written as X[o,a] ( c (t — £2) + Trd(k — n)/n). The combi- 
nation of these two representations with a logical product 
(AND) yields the following representation: 

fk{t;h,h) =X[o,a\ f c (* ~ 

■X[o,a] (c(t-t 2 ) + ^(fc-n)) . (3) 

The function /fe(i;ti,i2) corresponds to the logical value 
of the output Xk- fk(t;ti,t2) = 1 means that the pulses 
collide with each other near gate Xk, and fk(t; ti,t 2 ) = 
means that the collision does not happen near gate Xk- 
In the representation in Eq. we have done a reduction 
and neglected the pair annihilation of the pulses. Note 
that Eq. |2| is a logical expression that explicitly includes 
the continuous variable of time t. 

As we stated above, the "direction detector" in Fig. 
n uses a time difference between two signals with no 
discretization of time by a ruling clock. With the use 
of such a "time-logic operation," it should be possible 
to construct a novel type of parallel computing ma- 
chine in the absence of a ruling clock. For a Neumann- 
type computer equipped with a CPU to deal with such 
"time-dependent" signals, time information has to be dis- 
cretized, recorded on a one-dimensional matrix and pro- 
cessed in a step- wise manner. On the other hand, living 
organisms continuously process many signals without any 
governing clock, and signals with time-sequential infor- 
mation are encoded in wave packets or successive pulses. 
It was reported that a type of owl uses an interaural 
time difference (ITD) to detect the direction of the sound 



made by its prey j2^,|2^|. Therefore, the time operation 
on an excitable field may more closely resemble the infor- 
mation processing in biological systems than Neumann- 
type computation. 

As we have described above, the promising properties 
of "field computation" can be realized in a chemical re- 
action system, i.e., the BZ reaction. However, the prob- 
lem with the BZ reaction is its fragility; the liquid reac- 
tion medium is not suitable for use in a practical com- 
puting machine. Thus, "field computation" that takes 
advantage of time-sequential information could be im- 
plemented in a solid circuit by using modern electronic 
technology, as in semiconductor production |3^| or in a 
solid surface [34| . To overcome this fragility, it may be 
necessary to expand the essential idea of "field compu- 
tation." Recent advances in computer development have 
involved the fabrication of networks, whose characteris- 
tics can be described by coupled differential equations or 
finite-difference equations without an explicit inclusion of 
real space. However, due to the many effects of miniatur- 
ization, it may be necessary to explicitly consider spatial 
extension at the individual nodes in a network. In other 
words, we have to describe the electronic network with 
partial differential equations instead of coupled ordinary 
differential equations. Thus, the novel idea of "field com- 
putation" may offer a path for the future development of 
"computers," by avoiding a serious bottleneck in current 
Neumann-type computation. 
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